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Abstract. This is a study of the scattering and absorption of planar gravitational waves by 
a Kerr black hole in vacuum. We apply the partial wave method to compute cross sections 
for the special case of radiation incident along the rotation axis. A catalogue of numerically- 
accurate cross sections is presented, for a range of incident wavelengths Adu) < 4 and rotation 
rates a < 0.999M. Three effects are studied in detail: polarization, helicity-reversal and 
glory scattering. First, a new approximation to the polarization in the long-wavelength limit 
is derived. We show that black hole rotation distinguishes between co- and counter-rotating 
wave helicities, leading to a term in the cross section proportional to aiv. Second, we confirm 
that helicity is not conserved by the scattering process, and show that superradiance amplifies 
the effect. For certain wavelengths, the back-scattered flux is enhanced by as much as 35 
times for a rapidly- rotating hole (e.g. for a = 0.999M at Moj = 0.945). Third, we observe 
regular glory and spiral scattering pealcs in the numerically-determined cross sections. We 
show that the angular width and intensity of the peaks may be estimated via a semi-classical 
approximation. We conclude with a discussion of the observable implications of our results. 
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1. Introduction 

Gravitational waves are propagating ripples in spacetime that are predicted by General 
Relativity (GR). There is strong indirect evidence for their existence, for example, from thirty- 
five years of pulsar timing measurements [1]. Due to the tiny amplitude of waves reaching 
Earth (with a dimensionless strain h ~ 10~^^), gravitational waves are yet to be detected 
directly. Many experimentalists are optimistic that first-light observations will occur within 
the decade. Ground-based detectors such as LIGO [5] are reaching sufficient sensitivities such 
that even non-detection will constrain astrophysical theories (for example, on the origin of 
gamma ray bursts [3]). Plans for a space-based detector (LISA) are well advanced jl]. 

Gravitational waves are of interest to astronomers because they are generated by some 
of the most energetic astrophysical processes, such as binary mergers, supernovae and galaxy 
collisions. Electromagnetic radiation carries relatively little information about the most 
energetic regions of such processes, because photons are strongly scattered, absorbed and 
thermalized by intervening matter. On the other hand, gravitational waves are only weakly 
coupled to matter, and carry information about the dynamics at the heart of such processes. 

Black holes are likely to be powerful sources of gravitational waves. For example, at the 
end of the inspiral phase, the merger of a pair of black holes creates a distinctive "chirp" signal. 
Rapid progress has been made in numerically modelling the merger process [5], providing 
template signals for experimentalists seeking to separate a weak gravitational-wave signal 
from a noisy background. 

In addition to emitting radiation, black holes will scatter and absorb radiation that 
impinges upon them. To some extent, this is a secondary effect. However, scattering processes 
have long been of foundational interest to physicists, and this process is no exception. In this 
paper, we seek to answer the following question: how are gravitational plane waves scattered 
and absorbed by rotating (Kerr) black holes? In our scenario, we assume a monochromatic 
wave produced by some far-distant source impinges on an isolated black hole. We assume the 
incident wave is essentially planar, long-lasting, and weak enough to be treated through the 
linearized approximation. Here, we consider only the special case of on-axis incidence. Our 
principle aim is to determine the absorption cross section o"a and the differential scattering 
cross section da/dQ as a function of scattering angle 6. Inevitably, the long-range 1/r 
"Newtonian" component of the potential leads to a divergent scattering cross section in the 
direction of incidence. In this respect, a black hole is no different to any generic massive source. 
More interestingly though, the signal scattered through large angles carries information about 
the near-horizon structure of the hole. In this paper we compute this signal and examine in 
detail the effect of black hole rotation. 

Many studies of black hole scattering have been made over the last forty years. The most 
relevant for our purposes are studies of time-independent scattering, where it is assumed that 
the incident wave is monochromatic and long-lasting. Pioneering contributions include the 
works of Hildreth [6], Matzner and co-workers [3 El El [TDl E], Sanchez [121 UHl El [15] 
and Mashhoon [161 EZl [H]- Early studies focussed on the simplest setup: scalar (Klein- 
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Gordon) waves incident on a non-rotating hole. Later works extend the analysis to 
higher spin: fermionic [19], electromagnetic |20] and gravitational [21] fields. A number 
of authors have considered the polarizing effect of a rotating hole on waves with spin 
[T6l [T71 [T8l l22l l23l EH [251 ES]. The physical origin of regular oscillations in cross-section- 
versus-scattering-angle plots was clarified by DeWitt-Morette [27], Zhang |28j and others 
[291 130] . Such oscillations are a diffraction effect due to interference of null rays passing close 
to the unstable photon orbit. In other words, the large-angle scattering cross section contains 
information about the near-horizon structure of the hole. 

The efforts of Matzner and co-workers [H IH HOl HH El] culminated in the publication 
of a monograph [31], which set out in detail the partial- wave approach to computing black 
hole scattering cross sections. In this work, a range of plots of numerically-determined cross 
sections were presented. We believe that the plots for the rotating black hole are not accurate, 
due to numerical difficulties arising from the long-range nature of the gravitational field. Our 
belief is supported by more recent work by Glampedakis and Andersson [32] and [33], who 
examined the scattering of scalar waves by a Kerr hole. A primary aim of this paper is 
to combine the numerical accuracy of Glampedakis and Andersson with the partial wave 
methods of Matzner et al. for the gravitational wave case. 

A further objective is to examine in detail the coupling between the rotation of the hole 
and the helicity of the incident wave. The black hole is expected to have a polarizing effect [16] . 
A key result of this paper is that the polarizing effect is present even in the long-wavelength 
limit. A further observation is that parity-dependence of the scattering potential gives rise 
to a non-zero amplitude for helicity reversal. The helicity- reversed flux is maximal in the 
antipodal direction. We find that, close to the extremal rotation limit, the fiux scattered in 
the backward direction can be significantly enhanced by the effects of superradiance. 

This paper is structured as follows. In Section |2] we outline basic results for black hole 
wave scattering and quote the key analytic results of this work. In Section |3] we introduce 
the partial-wave method of Matzner et al. ( |3.1[ ), Teukolsky's equations (3.2), and the spin- 
weighted spheroidal harmonics (3.3). In Section [i] we apply the partial wave method to derive 
a new approximation for the scattering cross section and polarization in the long wavelength 
regime. In the process, we derive asymptotic results for the phase shifts (4.1), the spheroidal 
harmonics (4.2) and the scattering amplitudes (4.3). Section |5] describes methods for the 
numerical calculation of cross sections. The numerical task splits into three parts. First (5.1) 
we compute the spin-weighted spheroidal harmonics using a spectral decomposition method. 
Second (5.2) we calculate phase shifts using a Sasaki-Nakamura transformation. Third (5.3) 
we introduce a series reduction method for improving the convergence properties of partial 
wave series. In Section |6] we present numerical results for the absorption and scattering cross 
sections and the polarization. We conclude in Section [7] with a discussion of the physical 
implications of our results. Note that units G = c = 1 are employed throughout this paper. 
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2. Scattering by Black Holes 

The scattering scenario is characterised by two dimensionless parameters, 

M\uj\ = nrs/\, and a* = a/M = J/M"^ (1) 

The coupling Muj expresses the ratio of Schwarzschild horizon rs to incident wavelength 
A, and < a* < 1 is a measure of the rotation rate of hole. Here, M and J are the 
mass and angular momentum of the hole, and uj is the angular frequency of the incident 
wave. Note that the frequency uj may take either sign, depending on the polarization of 
the incident wave. Positive > corresponds to a circularly-polarized incident wave co- 
rotating with the black hole, whereas negative < corresponds to a counter-rotating 
circular polarization. Approximations for the scattering cross section are available in two 
limiting regimes, M\uj\ ^ 1 and M\uj\ > 1. 



2.1. The Long-Wavelength Regime, M\uj\ ^ 1 

If the wavelength is much larger than the horizon size, M\uj\ <^ 1, then the black hole 
interaction may be treated through perturbative methods. Various authors have shown 
[221 EH 133 ESI EZ] that, in the non-rotating case (a = 0), the limiting cross section depends 
on the spin s of the (massless) wave as follows. 



Mu^-.o dnj 



sin4(6»/2) " 
cos^(e/2) c - 1 /9 

sin4(e/2) ^ , . 

sm4(6»/2) 

where 6 is the scattering angle. Regardless of spin, all cross sections diverge in the same 
manner in the forward direction, da/dQ ~ 16M^/^'^. The divergence is due to the long-range 
nature of the gravitational potential. Through semiclassical arguments, it may be argued that 
the divergence is a consequence of the weak-field Einstein deflection (i. e. in the far-field, an 
incoming geodesic with an impact parameter b is bent through an angle 6 ~ 4M/b). 

It is interesting that all cross sections in ^ obey the general rule da/dVt ~ 
cosl^'^l(6'/2)/sin'^(^/2), except for the gravitational wave cross section [351 EZj; which 
includes an extra term. The extra term is due to the parity-dependence of the interaction of 
the linearised wave with the black hole. Odd-parity (axial or Regge- Wheeler) and even-parity 
(polar or Zerilli) partial waves with the same I pick up different phase shifts on scattering. 
This leads to the non-conservation of helicity in the scattering of a (linearized) gravitational 
wave. 

If the black hole is rotating, then an unpolarized incident wave will acquire a partial 
polarization. It is shown in Section |4] that the cross section for a long- wavelength gravitational 
wave impinging along the axis of incidence of a rotating hole is 

^da cos^(9/2) r , sin®(6'/2) r o,^, ,i 
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The cross section (|3| depends on the sign of u, so scattering from a rotating black hole induces 
a partial polarization, 

^ = f(o;>0) + ^(.<0) - ^^/'^ lcos«(^/2) + sin«(^^/2) J " 

In the small-angle limit, the polarization V ~ —\uj\{J/M)6^ is half that predicted by 
Guadagnini and Barbieri [231 [2H [2S1 12S] for scattering from classical rotating matter. 

2.2. Glory and Spiral Scattering, M\uj\ > 1 

When M\uj\ > 1, interference effects arise in the scattering cross sections [29]. Rays which pass 
in opposite senses around the black hole interfere to create constructive or destructive fringes. 
In the backward direction, a 'glory' will arise. This is a phenomenon familiar from optical 
physics [3H]- Strictly speaking, a glory is a bright spot in the backward-scattering direction 
{6 = it). However, we will use the term more loosely to describe the interference pattern near 
9 n. Using semi-classical arguments, DeWitt-Morette, Zhang and others [271 [2H1 |29l [30] 
derived an approximation for the scattering cross section close to the backward direction, 

db 

— J2\s\'^{bgUJsm9), (5) 

9=TT 

where b is the impact parameter, and bg is the impact parameter for which the deflection 
angle is = vr. Here J2\s\ is a Bessel function and s is the spin of the field. It is worth 
noting that ([s]) is zero in the backward direction if the perturbing field has spin (hence, a 
bright spot on-axis is absent). The gravitational- wave glory appears as a ring around 6 = n, 
with an angular radius proportional to 1/M|ci;|. In addition to a glory ring, we expect to 
see "spiral scattering" oscillations at intermediate angles [30,J. These oscillations arise from 
the interference between overlapping classical paths that are scattered through the angles 6, 
27C -e,27T + e, etc. 

To extend the glory-scattering approximation (jSj) to the rotating hole it is necessary to 
compute the glory impact parameter bg and the derivative db/dO for scattering geodesies. In 
the Schwarzschild case, these quantities may be computed by numerically solving the orbit 
equations, or via the deflection-angle approximation of Darwin [39] , 

-.nf^V (6) 



^3.48M^ 

This approximation is valid for impact parameters b sufficiently close the critical value b = be. 
We find the exact values are bg = 5.3570M and b'^\db/d6\ = 4.896M^ whereas the Darwin 
approximation gives bg = 5.3465M and b'^\db/d9\ = 4.30M'^. 

Let us now consider incoming geodesies that are initially parallel to the rotation axis of 
the hole. The defiection angle 6 is determined by the integral equation [40] 

(7) 



Joo VR 
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a* = 


0.0 


0.25 


0.5 


0.7 


0.9 


0.99 


1.0 


rjM 


3.00000 


2.97190 


2.88322 


2.75791 


2.56000 


2.43100 


2.41421 


bc/M 


5.19615 


5.17791 


5.12053 


5.04015 


4.91606 


4.83828 


4.82843 


C{a) 


1.000 


1.010 


1.042 


1.094 


1.200 


1.293 


1.307 


D{a) 


3.482 


3.461 


3.384 


3.247 


2.942 


2.665 


2.623 



Table 1. Numerical values of the unstable orbit radius (re), the critical impact parameter 
(6c) and the coefhcients in the logarithmic approximation to the scattering angle (C(a) and 
Dia)). 



where 



R 

e 



r'^ + ^2a'-by^ + 2b^r-a\b'-a^] 
b^-a^sin^e. 



(8) 
(9) 



For b > be, the quartic R has four distinct real roots Tq < ri < r2 < r^, of which (the radius 
of closest approach) is the largest. For the critical case, b 
at the photon orbit radius (rc = r2 

1 



be, the roots r2 and coincide 



— 



((27M^ 



r3 and dR/dr\r=r^ = 0). Explicitly, 



Te = v/3 + (3M2 - a^)/v + M, 



(10) 

(11) 
(12) 



= 27M{M^ 

In Appendix A we derive a logarithmic approximation similar to ^ for the polar orbits that 
pass close to the critical radius. We find that 

b-be 



9 ^ -C{a)ln 



(13) 



D{a)M^ 

where C{a) and D{a) are dimensionless constants which depend on the rotation rate of 



the hole. Closed-form expressions (A. 10) and (A. 11) for C{a) and D{a) are somewhat 
complicated; for convenience, numerical values are listed in Table [T] for a range of a^,. 

We may use the logarithmic approximation (13) to estimate bg and dbg/dO in Eq. (|5|. 
For = 0.99, the logarithmic approximation gives bg = 5.0609M, dbg/dO = — 0.1722M, 
and b'^\dbg/d6\ = 4.4109M^ which may be compared with exact values bg = 5.0925M, 
dbg/de = -0.2209M and bl\dbg/de\ = 5.728QM^, obtained by solving Q numerically. Hence 
(13) is a reasonable but not precise approximation for our purposes. In Fig. 17 the semi- 



classical approximation is compared with exact cross sections calculated via the partial wave 
approach. 



3. The Partial Wave Method 



To compute absorption and scattering cross sections, we invoke the partial wave methods of 
Matzner et al. [51 [IDl HB EI] • Since the full analysis is somewhat involved and is presented 
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in full elsewhere we just recap the key results here. We restrict attention to the special 
case of a gravitational wave impinging along the rotation axis of a black hole. 

3.1. Partial Wave Series 

The scattering cross section is found from the sum of the square of two amplitudes, 
The amplitudes may be expressed as partial wave series, 



P=±l 1=2 

oo 

9(0) = - E E ^(-1)' -2^'(0' -2^'(^ - ^' [^^P (2<J - 1] • (16) 



P=±l i=2 

In these expressions, exp(2i5;'^^) are phase factors which must be determined from a radial 
equation, _2Sf{9]auj) are spin-weighted spheroidal harmonics [H], and 9 is the scattering 
angle. Note the presence of a sum over parities, P = ±1. 

In the far-field limit, the metric perturbation /i^iT'* expressed in the outgoing transverse- 
traceless gauge takes the simple asymptotic form |H] 

h^^J'^m^'m'' ~ ^ — + , (17) 



where r* is the radial tortoise coordinate defined in (27), and m'^ is a (complex) null vector 
lying on the surface of a unit sphere. f{6) represents a helicity-conserving amplitude, whereas 
g{6) represents a helicity-reversing amplitude. The helicity-reversing amplitude is non-zero 
in general, because the phase shifts depend on parity P. 

The absorption cross section cTq may also be computed from the phase shifts. It is given 
by a sum over angular modes, 

/I 2 °° 

cTa = ^y^\-2Sf{0;auj)\'Ti2, (18) 

1=2 

where < Tim < 1 is the transmission factor for a given mode, 

2 



Im 



(19) 



Hence absorption is related to the imaginary component of the phase shift 6f 



3.2. Teukolsky Equations 

In the early 1970s, Teukolsky and Press |l2l SSj HH Il5] applied the Newman-Penrose formalism 
P6] to derive perturbation equations for the Kerr black hole. Somewhat unexpectedly, they 
found that a complete separation of variables was possible for the Weyl scalars ^'j. For 
instance, 

p-4v|>4 = e— * e™<^ _2Sr{0; au) _2i?/™.(r). (20) 
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With this ansatz, the Newman-Penrose equations decouple into radial and angular ODEs, 

^ Tr ^ * + k ^+*...-A,J fl = (21) 



and 



1 d f „dS\ /oo 9^ m'^ 2ms cos 9 



sinOde V de ) ' V sin^^ sin^ ^ 

-2auj s cos e - s'^cot^e + s + AiJj S = Q (22) 



where 



A = - 2Mr + (23) 

K = (r^ + a^)uj - am, (24) 

km = Am + d^^'^ - 2amu), (25) 
and s = —2 in this case. 



The angular equation (22), together with boundary conditions of regularity at = and 



6 = Tx constitutes a Sturm-Liouville eigenvalue problem. Hence the eigenfunctions — the 



spin-weighted spheroidal functions that appear in (15) and (16) — form a complete basis. A 
method for their calculation is outlined in Section 13. 3[ 

The Kerr black hole has two horizons, at r± = M±\^]\/P — a^, where A = 0. The s = —2 
radial equation is solved subject to the condition that the radiation flux is purely ingoing at 
the outer horizon r+, 

lim _2-R/mc<; ~ A^ exp(— z(co' — mr2/i)r*), where fi/i = a/(2Mr+), (26) 

r— >r-f 

and r^, is the tortoise coordinate defined by 

dr^ _ (r^ + a^^ 
dr ~ A 

iti 

lowest order 



(27) 

Towards spatial infinity, there are two linearly-independent solutions, -R;™^^ and ■ To 



Rtt = r-' expi-tu;r^) [1 + 0{r-')] 
Rt^u = exp(z^r,) [1 + 0{r~')] 
Hence the ingoing solution obeys the boundary conditions 



^(trans) s p~i(uj~mQ.fi')r 



r* — OO 



where B'fl^^\ and B^^ are complex constants. The phase shifts appearing in (15) 



and ( 16 ) are defined by |TOl 

exp(2,0-(-ir( "-'"\;:r-^ )^ (30) 



Imu} 
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where 



[Re(C)]^ 



((A + 2)2 + Aamu - Aa'^u'^) [\^ + 36amcj - SGa^cu^] 
+ (2A + 3) (960^0^2 - ASaujm) - lUuj^a'^. 



(31) 



Note the exphcit parity-dependence in (30). 



A method for finding analytic solutions to Teukolsky's radial equation (21) in the long- 



wavelength regime {Mu <^ 1) is described in section |4| Outside this regime, we will use 
numerical methods. In principle, to compute phase shifts all we need to do is integrate the 
radial equation out from near the horizon to large r, where we may read off the coefficients 



B;Z! and 5, 

the factor of difference between ingoing and outgoing asymptotics — causes great numerical 
difficulties. A resolution to this problem is given in Section [5j 



(refl) 



'imu! "''-'■'^ ^imJ ■ -'■^ practice, it is not so straightfoward, since the 'peehng' behaviour [16] 



3.3. Spheroidal Harmonics 



A range of methods [41j exist for solving the angular equation (22). For instance, to compute 



the eigenvalues \im we might employ either small-ati; [47j or large-aa; [IE] approximations, 
a continued-fraction method |19], or a finite difference method [50]. Here we make use of a 
spectral decomposition approach described in Appendix A of [51]. 

Following the spectral decomposition approach, each spheroidal harmonic is expressed as 
a sum of spherical harmonics. 



,Sr{0] auo) 



E 



.(0 



(32) 



j=max(|?n|,|s|) 



where are expansion coefficients to be determined. Note we have suppressed the azimuthal 
dependence of the spherical harmonics s^™'{9). The spherical harmonics satisfy the equation 



sin 9 dO 



d ( . d,Yp 



+ 



J(J + 1)- 



+ 2ms cos ^ + 
sin^ 9 



0. 



Hence, by substituting ansatz (32) into the angular equation (22) one gets 



(33) 
(34) 



J2 cos^ 9 - 2aws cos 9 - j{j + 1) + Ei^] ,F™(^) = 0. 

j=max(|m|,|s|) 

where Ei^ = Aim + s{s + 1). Now, we take advantage of the orthonormality of the spherical 



harmonics, / d^l gY,^ (9) (9) = Skj. Multiplying (34) by sY'^™'(6') and integrating yields 



2 2 (2) 



+ 



+ 



2 2 (2) 



2aujsc 



(1) 



dl) 



+ 



(1) 



2 2 (2) 



2aujsc 



(1) 

k,k+l 



^k+1 



+ 



2 2 (2) 



^k+2 



-Elmbf. 



(0 



where 

.(2) _ 



-kj 



dnsY^{9)cos'^9sYl^i9) 



dn sY^ (9) cos 9 sY"" (9) 



^kj 



2j + l 
3 V 2A; + 1 



(j, 2, m, 0\k, m) {j, 2, -s,0\k, -s) 



2j + l 
2A; + 1 



(j, 1, m, 0\k, m) {j, 1, -s, 0\k, -s) . 



(35) 

(36) 
(37) 
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The numbers {ji, j2,^i,fn2\j,^) are Clebsch-Gordan coefficients [SB]. 

Equation (35) can be rewritten as a matrix equation, with the numbers bf' representing 



the coefficients of the matrix's eigenvectors, and —{Aim + s{s + 1)) representing the matrix's 
eigenvalues. The matrix is band-diagonal, so the task of finding the eigenvalues and 
eigenvectors is straightforward and numerically efficient. 



To calculate the spheroidal harmonics from (32), we need a method for computing the 



spherical harmonics of spin-weight s = —2. A suitable recursive method is described in 



Appendix B 



4. Polarization of Long- Wavelength Waves 

The infiuence of rotation on the scattering of long- wavelength waves (M|a;| -C 1) has been 
considered by a number of authors P [161 113 HI ISSl [231 1211 [251 [26]. It seems that no clear 
consensus has yet emerged. Considering electromagnetic radiation, Mashhoon [TB] noted that 
"one must expect partial polarization in the scattered light when an unpolarized wave is 
incident on a Kerr black hole", and went on to observe that "the polarizing property of a 
Kerr black hole is probably maintained for very low frequencies" . However, the corollary that 
incoming gravitational waves will also be polarized at low frequencies was thrown into doubt 
by a study [22] in which the scattering amplitude was computed using Feynman-diagram 
techniques. The authors concluded that "the angular momentum of the scatterer has no 
polarizing effect on incident, unpolarized gravitational waves" (even though they found that 
unpolarized electromagnetic waves were polarized). Some doubts about the gauge-invariance 
of the results in have been raised. To clarify the issue, the Feynman-diagram approach 
was recently revisited and improved in p6]. In a series of papers, Guadagnini and Barbieri 
[231 IMl [251 126] have argued that the net polarisation V induced by rotating classical matter 
is 



^(^>0) + ^(c.<0) 



T ^ t ) ^ n , f ■ .r!. ^ -\sHJ/M)e' (38) 



where 6 is the scattering angle, which is assumed to be small. 

In this section, we apply the partial wave method described in the preceding section to 
derive a long-wavelength approximation to the polarization of gravitational waves impinging 
along the axis of a Kerr black hole. The key results of this analysis were presented in Eq. ^ 
and Curiously, we find the polarization at small angles to be exactly half that given by 



(38). 



Our objective is to compute the scattering amplitudes / (15) and g (16) to lowest order 
in Mu) and auj. In [37] it was shown that, to lowest order in Mu, the Schwarzschild (a = 0) 
amplitudes may be written as 

/schw(x)-Me r(l + ,e) [|(l-x)]i- ^^^^ 
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and 

gscUx) = - ^)] (40) 

where $ = 2eln |2e|, x = cos 6, e = 2Mu and T{y) is the Gamma function. In the following 
analysis, we show that black hole rotation (a* > 0) introduces a new term in each amplitude 
which is proportional to au. 



4.1. Phase shifts via the Mano- Suzuki- Takasugi method 



Low-frequency analytic results for the phase shifts (30) may be obtained from the formalism 
developed by Mano, Suzuki and Takasugi [Ml [55] (MST). In the MST approach, reviewed 
in [56], solutions to the Teukolsky radial equation satisfying boundary conditions (29) are 
expressed as infinite series of special functions. Two series are used. The 'horizon' series of 2-F1 
hypergeometric functions is convergent at all radii up to (but not including) spatial infinity. 
The 'far-field' series of Coulomb wavefunctions is convergent up to (but not including) the 



outer horizon. The complex constants B 
by matching the two series. It is found that 



(inc) T5(refi) 



(trans) 



^(inc) 

I 

B 



(refl) 



Imui 



Kjs) - ie~ 



,Sm TT z/ 



defined in ( 29 ) are determined 

s + ie)) 



(41) 



-l-2s 



sin (7r(z/ + s — ie)) 

^00 ^ - [K,{s)+te''"'K^,.i{s)] A''_. (42) 

Here Ag is just a normalization constant, and we follow the conventions of [SI] by defining 

e = 2Muj and k = ^1 - aj. (43) 

Note that e must be positive in the MST expressions; the u; < results may be found 
via the symmetry of the radial function sRimw = sR*i-m~w from which it follows that 

p5{inc/rcfl) Tj(inc/refi)* 

lm~u) l—mu> 

The parameter u is known as the "renormalized angular momentum" and has the low- 
frequency expansion 



2/ + 1 



+ 



21 2 



{P - s 



2\2 



+ (2/ + l)(2/ + 3)(2/ + 3) (2/ - l)(2/)(2/ 

Note the absence of a linear term in e. 

The coefficients K^y and K_i,_i may be computed via a complicated series expansion, 
detailed in [5l]. This is not necessary for our purposes, because K^y^ilK^ ~ (9(e^'~^). Hence 
K_y_i may be neglected, and the terms will cancel upon taking the ratio of b\^^ and 

p(refl) 



The coefficients AK and A'^_ are given by 



A^ 



2-l+Sg-iegi(7r/2)(i/+l-s)g-7r£/2 



It 



Tiv + l 



le 



E 



eg-i{7r/2){i/+l+s)g-7re/2 



E 



[v + l + s 



le 



[v + 1 — s + ie), 



(45) 
(46) 



Scattering and absorption of gravitational plane waves by rotating black holes 



12 



where a'^ ~ (9(e'"''). Here, we wish to conduct an expansion accurate to second order in Mco; 
hence we require exphcit formulae for a'^2y ■ ■ ■ ,0'2- These were computed in [51] and are hsted 
Appendix C 



m 



After inserting equations (C.l - C.5 ) into (30 ), and taking some care with series expansions 

r(/ 



in the small parameter e, we find 



g2ie In 2e ^—leK ^—i-K(v—l) " 



1 



1 + + 



(47) 



for positive uj. Note that there is no first-order term in the square brackets; the linear terms 
cancel exactly. The second-order coefficient is 



a 



(2) 



12ima^ 



(48) 



'"^ /(/ + 1) (/- l)/2(/ + l)2(/ + 2)' 

where m is the azimuthal number. Note that the e"*'^*^'""') and e^*^/'*^'"*"^) factors also give an 
/-dependent contribution at second order in e^. 



Via the symmetry -Sf^'^™^'' 



)(ref/inc) 

Hence the phase shift may be written more generally as 



^i-mil'^'" follows that exp{2i6^ 



Im—uj' 



exp(2z5j" 



g2je In |2e| g~*^'^g~* ■k{v—1)u)/\uj\ 



T{l + l 



le 



I + a'S^e' + 0{e')] . (49) 



This expression is valid for either sign of e = 2Muj. 

In a previous study [57|, the low- Ma; approximations of Poisson and Sasaki [SZ] were 
used to show that the Schwarzschild (a* = 0) phase shifts are 



exp(2i5, 



IrriLu) 



g2ie In 2t^—ie 



r(/ + 1 + ie 



(50) 



This is consistent with (47) since k = 1 in the non-rotating case. Moreover, it is remarkable 



that, when a^, = 0, equation (50) actually holds to one order higher in e than first supposed. 



4-2. Series Expansion of Spheroidal Harmonics 



In order to compute low-frequency scattering amplitudes (15, 16) we must first expand the 



spheroidal harmonics to second order in the spheroidicity parameter z = auj. Noting that 
&;2n ~ C'(z"'), let us make the expansion 



"l-2 



2^(0) 



2^(1) 



z'd 



^2) 



"+1 
_ 2 .(0) 



The normalisation condition implies that (c/q) 

7(0) ^(0) ^(1) ^(2) ^(0) ^(1) ^(0) 



+ «i) 



(0)x2 



(51) 
(52) 
(53) 

The remaining six 



unknowns {c?_2, (r_{^(r_[^ (Xq , (il^i, (i^_2} determined from the equations 



(/ + \)d\ 
2{l + l)d^_ll 



-so 



(1) 



SC, 



(1) 



'-2sc 



(1) 

1+1,1+1 



-^/mJ"+l + (^1+1,1 



(54) 
(55) 
(56) 
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(1) 
-1 



(2) _0_(1) ^(0) 



^(1)^^(0) 



.(2) 



2(2/ + 3)d^°^ 



+2 



-2(2/ - l)d 



(0) 



'^1-2,1 ^■'''/-2J-l"-l' 



where the angular eigenvalue has been expanded in powers of z = aut as 



fc=l 



and the first few coefficients are determined by the identities 

^(1) 



(2) 

n 



(3) , J2) ^(0) , (2) ,0 _2sc(i) (7«-2sc^^^ rf^^^ 



Im 
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(57) 
(58) 
(59) 

(60) 



(61) 
(62) 
(63) 



For this calculation, we need only the m 



-2 harmonics. The required 



Appendix C 


.( 


C.6 


C.13 



_2, d';^{, d]^l, d^Q \ (i^l, c/^], c/1^2}) which are listed in 
To compute the spheroidal harmonics explicit 



Gordon coefficients into (54 59|) yields explicit ex pressions for t he expansi on coefficients 
{^(0) ^(0) ^(1) ^(2) ^(0) ^(li^ftr 



Appendix C, (C.14-C.19). 



y we require expressions for the spin- 
cos^. These may be found by acting 



weighted spherical harmonics _2l^^(a;), where x 

on spherical harmonics of spin- weight zero, qY^[x) = Pi{x), with ladder operators 

The spin- weight is lowered with the operator 6, and the azimuthal number is raised with L+. 
These operators are defined by 

m + sx 



s.Yrix) 



Vl — X2 (9^ 



x" 



,Yr{x) = -^{i + s){i-s+i),_,Yr{x), 



(64) 



s + mx 



x^ 



sYi^'ix) = ^/{l-m){l + m + l)sY{^+\x). (65) 



Here, is shorthand for the partial derivative with respect to x = cos^^. By acting with 
(5L+(5L+ on oYi^{x), it is straightforward to show that the spin-weighted harmonics in (15) 

2/ + 1 LPi{x) 

' V 47r (/ - 1 



and ( 16 ) can be written 



lY.'ix) 



An (/- l)/(/ + !)(/-!- 2)' 

'^^^ x)d,d^,{l-x)d,Pi{x). 



[l + xVdJl 



(66) 
(67) 



Their values in the forward and backward directions are particularly simple, 

-1) 



-2Ynx=l) 



21 + 1 



4tt 



YHx 



-2^1 



0. 



This implies that the values of the spheroidal harmonics in the forward direction are 

.,sr{x = 1) = 



2/+ 1 



4tt 



(68) 



(69) 
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5o 



^ 16 ^ 192(/2 + / + l)^ 
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(70) 
(71) 



(2/ + 3)2(2/- 1)2 + (/ + 1)4/4 (/ + 1)4/4- 

4-3. Scattering Amplitudes 

Let us now calculate the helicity-conserving and helicity-reversing amplitudes (/ and g) in 
the long- wavelength limit. 



4-3.1. Helicity-conserving amplitude. First, consider the helicity-conserving amplitude / 



defined in (|15|). To begin, note that the sum of the positive and negative-parity phase terms 

r(/-i-k)r(/ + 3) 



can be written 



P Imuj -I- P lmu> 



3«X 



2 V' ' ' J r(/ + 3 + ie)r(/-l) 

where e*^ = e2iein|2e|g-ieK second-order coefficient is 



1 + P?2^' + o{€^ 



,(2) 
Im 



UJ f V — I 



\uo\ 



%— + 
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/(/ + 1) /(/ + 1) (2/ + 3)(2/- 1)' 



(72) 



(73) 



Inserting the expansions of the spheroidal harmonics (|32, 69), and using (72), we may write 
the / amplitude as 



/(^) 



3«X 



2iij) 



(74) 



where 
F{x) 

and 



2/ -3 (/ + l)(/ + 2) 
2/ + 1 (/-3)(/-2) 



;2/ + 3(/-l)(n 2/ + 5 {1-1)1 (0 

+ 27Ti(7+3)^'+i^'+^^'^^ + V 27n(I+W+4) '^'^ ^ 



(76) 



and l>f^ are the expansion coefficients defined in (32). 



To compute the higher-order terms in the s um F{x), we substitute in the explicit forms 

Next, we rewrite the sum so all terms 



for the coefficients bf^ calculated in subsection 
have a common factor of Pi. To demonstrate this process, let us begin by considering just 
the linear term in tu, which we denote F(„,\. We find 



4.2 



1=2 



r(/ + 3 + i£) I V 2/ + 1 (/ - 2)' 
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+ O{uj 



2> 




It is straightforward to verify that the terms in parantheses on the second hne cancel exactly, 
and so a linear- in-cj term in F{x) is not present. 



Note that to obtain P/(x) as a common factor in (77) we redefined the summation variable 
(i.e. / ^ / ± 1). Care must be taken in this process, since it changes the lower limit of the 
sum. In the above, we added an 'extra' P2{x) term to the series without consequence, because 
the coefficient d^^l{l — 1) is zero for I = 2. We also neglected a Pi{x) term (present in the 
top line but not subsequently). This is justified for our purposes because, to obtain the final 
amplitude /, we act on F with C^, and C^Piix) = 0. 

Let us now repeat this process and focus only on the part which is linear in a. This time, 
we will keep terms up to second order in u. Such terms arise from (i) the a^,-dependent part 
of the phase P^^e'^, and (ii) the linear term Si coupled to an e term arising from the effect 
of redefining the summation variable. We may split the result into two parts: a sum which 
turns out to be zero, and an / = 2 term which arises from the change of summation variable. 
That is, 

F{x) ^ Foix) + zt (Ff ) + F^=')) + 0{z\e^) + 0{J'). (78) 
The lowest-order term 

was defined and summed in [37j. The additional terms are 

1=2 



(/-!)(/ + 3) 
= + O(e), 

and 




-2(/ + l) /(2/ + 3) (0),, , ^, 4 , 21 (21-1) M), 



d'-'' (I + 1) - — + — J frfV (^ - 1 

' /(/ + !) (/-2)(/ + 2)V(2/ + l) 




Fr\-) = -2.^^P.(x). (81) 
r(5 + le) 

(1=2) 

The term zeF{ (x) is responsible for the lowest-order polarizing effect. It leads to a term 
in the scattering amplitude which depends on the sign of uj. That is, / ~ /schw + /poi where 

n ^2 

/poi = -2Maiue'^ 



1(1 + 



^(l + 0(.)). (82) 
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and /schw was defined in (39). Since /poi is in phase with /schw, this gives a first-order 

.cos«(^/2) 



contribution to the cross section, 

2|/schw/poi| = -4ao.M^ 



(83) 



4-3.2. Helicity-reversing amplitude. We now repeat the analysis for the hehcity-reversing 
amphtude, g{0)- To simphfy matters, we only expand the relevant terms to first-order to 
recover the polarizing correction. Let us begin by writing the amplitude as 



9(0) 



2iuj 



C^,G{x), where G{x) = Y,i-'^yWiVi{l + S^z + Oiz^)), (84) 



1=2 



and 

1^.^6.eS^^li^Mll 



Vi 



T{l + 3 + ie) r(/ + 3) 
Pi(-x) + z 



WTT] 



e + 0{e' 



^5) 



21 



^^ + ^rfL1(/)P.-l(-x)+.l/^^ + ^^~^'^(°)^-" ^ .^,....2, 



-d':>m+i{-x)+o{z').{86) 



2/ + 1/-2 ' ' ' y 21 + 11 + 3 
As before, we may move terms of the series up or down (/ ^ / ± 1) to obtain a common factor 
of Pi{—x). As before, we find that, to lowest-order, the terms in the sum cancel out, leaving 
only an / = 2 term which arises from redefining the summation variable. That is, 

G{x) ^ Goix) + ze f + G^^A + 0{z\ e^) + 0{uj^) 



and 



Gf\x) 



r(/-i-ie) r(/-i) 



1=2 



V{l + 3 + it) r(/ + 3) 



{21 + 1) 



16 



2(/ + 2)2 2(/-l) 



/(/ + 1) 



/2 



(^ + 1)^ 



(87) 



Pi(-x) 



G'i 



a=2), 



X] 



+ O(e), 
Til -it) , , 



^9) 



The zeG\ (x) term is responsible for the lowest-order polarizing effect. It leads to a term 
in the scattering amplitude which depends on the sign of a;. That is, g ~ (y'schw. + 5'poh where 

n 2 



X] 



r 1 



V{\ + %t) 



(90) 



and (7schw. was defined in (40 ). The first-order contribution to the helicity-reversal cross section 



IS 



2 l^schw^poil = 4awM2 sin^(^^/2). (91) 

Combining results (83) and (91) leads to the cross section ([3]), and hence the polarization (|4]). 
The consequences of polarization are discussed further in Section [7) In the next section we 
outline a numerical method for moving beyond the long-wavelength regime. 
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5. Numerical Method 



The task of computing scattering and absorption cross sections numerically may be divided 



into three steps. First, we must solve the angular equation (22 ) to determine the spin- weighted 
spheroidal harmonics _2Sf{0; au) and their eigenvalues Xim- Second, we must solve the radial 



equation (|21|) to determine the phase shifts e^*'^'^'^. Third, we must find a way of improving the 



convergence properties of the partial wave series (15 16) for a numerical calculation. These 



steps are described in sections 5A^, 5]2]and 5^ below 



5.1. Spheroidal Harmonics 

The the spheroidal harmonics -2Sf 
using the matrix method detailed in Section |3.3 



au) and angular separation constants -2Xim were found 
The accuracy of the eigenvectors of 



coefficients b^[^ was checked by back-substitution into (35). The spherical harmonics _2yj^(^) 
calculated via the method in [Appendix B| were verified to be highly accurate up to / ~ 100. 



The spheroidal harmonics -2Sf{6,auj) were tested against the results of Leaver's expansion 
|i9] . The eigenvalues were tested against Leaver's continued fraction method, and against 
low-aa; expansions jH]. 



5.2. Phase Shift Calculation 



As we saw in (28), the outgoing part of Teukolsky's radial function -2Rimuj dominates over the 
ingoing part by a factor of in the far field. This behaviour is due to the long-ranged nature 
of Teukolsky's equation. A solution to this problem was supplied by Sasaki and Nakamura 
[501 156] , who found a transformation that maps the Teukolsky function R{r) to a function 
X{r) governed by an equation with a short-ranged potential. 



To compute phase shifts (30), we take the following numerical approach. First, we find a 



power series representation for the Teukolsky function near the outer horizon. Then, starting 
at r = r+ + 7], where t] ~ 10~^M, we integrate the Teukolsky radial equation (21) out to 
r = Tm- Typically, ~ 5M. At Vm we transform to the Sasaki-Nakamura formalism and 
integrate (102) out into the far field, up to r = Too, where r^o ~ 



50M. There, we read off 
the ingoing and outgoing coefficients and thus compute the phase shifts. The phase shifts are 
insensitive to changes in rj, and Too. The details of this method are outUned below. 



5.2.1. Horizon expansion of the Teukolsky function. Teukolsky's radial equation (21) has 
a regular singular point at the outer horizon, r = r+. Using the method of Frobenius, the 
solution near this point may be expressed a power series in the distance from the (outer) 
horizon. Introducing the radial variable rj 



A{ri) dR B{r]) 
drf rj drj rj"^ 



R = 



r — r^, the radial equation (21 ) may be rewritten 

(92) 
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where R = -2Rimu- The functions Airj) and Birj) are 

-<^' - (-) 

_ + ^{r, + v)K - {8tu;{M + u + t]) + X) r]{7] + 2v) 

(r7 + 2z/)2 ' ^''^> 

where v = (r+ — r_)/2. Note that A{ri) and B{ri) are regular at = 0, so may be expanded 
in Maclaurin series, 

oo oo 

A{r^) = J2^kv', B{7^) = J2Bkv' (95) 

k=0 k=0 

The coefficients Aj^ are straightforward to compute: Aq = 1 and = (— l)'^"'"^/(2z/)'^ for 
k > 0. The coefficients 5^ are less simple, but are easily found with the help of a symbolic 
algebra package such as Maple. The radial function may be expressed as a power series. 



oo 



i?(r^) = ry7^i^,r7^ (96) 

k=0 

where 7 and Rk are constants to be determined. The index 7 satisfies an indicial equation, 
7(7 — 1) — 7^0 ~ Bq = 0, which has two solutions, 

2Mujr+ — am 

'J = za, and ■y = 2 — la, where a = . (97) 

r+ — r_ 

The latter choice corresponds to an ingoing wave. The coefficients Rj may be found from the 
recurrence relation 

j 

[(7 + j) (7 + J - 1) - (7 + j)Ao - Bo] R, = J2[{j + j- k)Ak + B,] R,.u. (98) 

fc=i 

Hence the series expansion may be taken to the desired order. 

5.2.2. Sasaki-Nakamura transformation. The Teukolsky function R{r) is transformed to a 
Sasaki-Nakamura function X{r) with the rule 

X{r) = [a{r)R + J (99) 

where 

air) = -iK{r)P/A'^ + 3idK/dr + A + 6A/r^ (100) 

P{r) =2A[-iK{r) +r - M -2A/r]. (101) 

The function X{r) satisfies the equation [50] . 
d'^X dX 

F{r)^-U{r)X = 0, (102) 



dr^ dr.. 



where the functions F(r) and U{r) are defined in Appendix D In the non-rotating limit 



a* = 0, the Sasaki-Nakamura reduces to the Regge- Wheeler radial equation j59| 156]. 
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Like the Teukolsky equation, the Sasaki-Nakamura equation admits two independent 
solutions at infinity, 



Xir) 



.{inc)p(inc)/ X 



I /l(refi) p(refi)/ X iccr* 



OO. 



(103) 

The principle advantage of the Sasaki-Nakamura formalism is that Pim^\r) and Pirn^\r) 
tend to unity as r oo. Hence it is straightforward to numerically determine the relative 
magnitude and phase of the ingoing and outgoing coefficents Af^^j and A\^^^ , which are 
(i-) D(-fl) in by 



related to 3^2! and B^^^ 



A 



(inc) 
Imu) 



A I ,2 R(inc) 



A 



(refl) 
Imu) 



_Co_^(refl) 



(104) 



where Cq is defined in (D.3) in Appendix D The phase shifts follow immediately from (30). 
The function Pi^^\r) has the series expansion 



A- B 



+ 



where |5T] 



(A + 2 + 2amu;) 



1 



B;^ = -- [(A + 2)2 - (A + 2)(2 - 4ama;) - 4 (ama; + 3iMa; - amuj{a7nuj + 2iMcc;))] 



6 



[Aamuj + i3in(A - 4 + 2amu + SzMcu) + 12M 



-2AnAMc<; - (ac<;)2(A - 3 + + 2amuo)\ . 

Whilst it has been stated (e.g. [51], eq. (4.21)) that P^^^\r) = P/i''^^*(r), this is not 
because F{r) and U{r) are complex. We find instead that 



•^out ^out 



out 



where 

+ CJCi/Co, 

-Bout = B*^ + [cj2c2 - u;ci(An + V2)] /co, 



+ 



[uorY 



In 



CJ"C3 — (^in + i)uj C2 + ( B\n + «An/2 ~ 2^ 2iMuj{aUJm — 1) ) UJCi 



/co, 



and Co, ci and C2 are defined in (D.3) in Appendix D This result is confirmed in 
erratum |52]. 



(105) 

(106) 
(107) 

(108) 
the case 

(109) 

(110) 
(111) 

(112) 
a recent 



5.3. Partial Wave Series Reduction 



The partial wave series (15) is divergent on-axis at 6' = 0. Physically, this is because the 
gravitational interaction is long-ranged; geodesies in the far-field are defiected through the 
Einstein angle 9 ~ AM/b. Since each term in the series is finite, an infinite number of terms is 
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required to reproduce the on-axis divergence. The magnitude of the series coefficients grows 
with /, and the series does not converge. 

Our computing power is finite, hence we must terminate the series at some / = 
/max- However, a naive truncation introduces unwanted 'noise' with a magnitude roughly 
proportional to 2Zniax + 1 and an angular width proportional to l//max- To avoid this problem, 
we employed a series reduction technique, inspired by a numerical method developed in the 
1950s [60] to compute Coulomb scattering series. 

Before applying the method, we first rewrite the partial wave series in terms of associated 
Legendre polynomials (and their derivatives). Series (15) is a sum over spheroidal harmonics; 
for brevity, let us write 



fi9) = J2fi-2Sf{6;au). 



;il3) 



1=2 



Using the results of Section 3.3, series (113) can be rewritten as a sum over spherical 
harmonics, 



f{9) = Y,Pj-2y'iO) where 

j=2 



(0 



1=2 



3.3 



(114) 
The spherical 



and b^j^ are coefficients found from the eigenproblem described in Sec. 
harmonics of spin-weight —2 may be expressed in terms of associated Legendre polynomials 
Pim{cos6) and their first and second derivatives. This decomposition is described in detail in 
Appendix B The series ( |113[ ) may be written 



/(^) = 

where x = cos 6 and 

JO) 



j=2 



cf\i-x': 



d'^p^. 



j2 



7(0) 



dx"^ ~^ ■' 



(0) 



dPj2 ^ ^(0) Pj2{x) 

dx ■' 1 — x'^ 



,(0) 



4(1 - x)AjFj, 



(115) 
(116) 

4vr (j-l)j(j + l)(j + 2)" ^^^^^ 
The idea behind the series reduction method is to improve the convergence properties of 
the series by removing the divergence at x = 1. We define a new series 



A, 



2j + l 



'l-xff\x) = Y,cf{l- 

j=2 



d^P. 



X 



j2 



dx'^ ~^ 



(k) dPj2 



dx ^ ^ 



(k)Pj2{x) 



I- x"' 



:ii8) 



ik) (k) (k) 

The functions Cj ,dj and are computed from (116) with the following recurrence 
relations. 



(k) 



d 



(k+i) 

'3 

(fc+1) 



(fc) 



(j - 2) 
2j-l 
(j - 2) 
2j-l 
(j - 2) 



Ik) 

S-1 



"7-1 



J + 3 (fc) 



2j + 3^+^ 

^4tW2(i+.)cr) 



jk) _ 

' 2j - 1 



(fc) 



J + 3 
2j + 3 



ef+\ + (1 + x)d 



(fc+1) 
i 



(119) 
(120) 
(121) 
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These relations were calculated with the aid of the recursive formulae (B.3, B.7, B.IO). In 



Appendix B we outline a stable recursive method for computing (1 — x^) '^^''•T , and 
The reduced series are found to converge more quickly as / — > cx). Typically, we apply 



two iterations of (119-121). In other words, we calculate the amplitude f{6) from (118) with 



k = 2. The series is terminated at /max ~ 60. 



6. Numerical Results 



In this section we present a range of numerically-determined absorption and scattering cross 
sections for a gravitational wave impinging along the axis of a rotating black hole. 

6.1. Absorption 

Figure [l] shows the absorption cross section aa (eq. 18) for a non-rotating black hole. The 
cross section is plotted as a function of gravitational coupling Mu. If the wavelength of the 
incoming gravitational wave is much larger than the horizon size (that is, Mu 1), then 
absorption is negligible. Conversely, if the wavelength is much smaller than the horizon (i.e. 
Muj ^ 1) then the cross section tends towards the geometric-optics limit of aa = 27ttM'^. In 
the intermediate regime Mu ~ 1, the contributions from successive angular modes (/ > 2) 
create a regular oscillatory pattern. 

Figure |2] compares Schwarzschild absorption cross sections for massless waves of spin 0, 
1/2, 1 and 2. The spin-0 and spin-1/2 cross sections approach o"a = IGttM^ and aa = 2ttM'^ 
as Muj [6T], whereas the s = 1 and s = 2 fields tend to zero in this limit. The cross 
section at low Muj is dominated by the lowest allowed angular modes (i.e. the / = |s| modes). 

Figure [3] shows the absorption cross section for an incident wave propagating parallel to 
the axis of a rotating (Kerr) black hole. In this plot, the incident wave is circularly polarized 
in the same sense as the rotation (i.e. > 0). In the low-Mcj regime, the cross section is 
negative. This is due to the superradiance effect which enhances the amplitude of the lowest 
corotating modes (particularly / = 2, m = 2). Hence, an ingoing wave may actually stimulate 
net emission from a rotating hole. 

In the high-coupling limit [Muj ^1), the cross section approaches the geometric-optics 
value, aa = vr6^. For a = 0.99M, the geometric cross section is approximately aa ~ 23.47rM^. 
The critical impact parameter be is found by solving ([T]) numerically; values for a range of a^: 
are listed in Table [H 

Figure |4] compares the absorption cross sections of co- and counter-rotating circularly 
polarized incident waves at a* = 0.99. It shows that a greater proportion of the counter- 
rotating helicity {u < 0) is absorbed than the co-rotating hehcity {u > 0). Superradiance 
is not stimulated by the counter-rotating helicity. In the limit Mu ^ 1, both cross sections 
approach the geometric-optics limit. 

Figure [5] illustrates how the absorption cross section depends on rotation rate a* for both 
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Figure 1. Gravitational wave absorption cross section of the Schwarzschild black hole. The 
top plot shows the gravitational absorption cross section aa as a function of gravitational 
coupling Mlli. The line at a a = 277rM^ shows the classical limit. The broken curves show 
the contributions ct; of each angular mode, from I — 2 onwards. The lower plot shows the 
relection factor for each angular mode. 



CO- and counter-rotating circular polarizations. For the former case, we find that the amount 
of stimulated emission created by superradiance increases as a* — > 1. Above Mu = 1, 
superradiance is prohibited, and the absorption cross section is strictly positive. For the 
prograde wave (cu > 0) it seems to be the case that aa{ai, Mu) < aa{a2, Mu) if oi > a2, for 
any Muj > 0. This is not the case for the retrograde polarization (a; < 0); the ordering of the 
curves is reversed as M\uj\ goes from zero to infinity. 



6.2. Scattering 

In this section we present plots of the differential scattering cross section da/dfl across a 
range of couplings Mu, and compare the non-rotating and rotating cases. 



6.2.1. Long wavelengths, Mtu <^ 1. In sections 2A_ and|4]we reviewed the analytic results 
available in the long- wavelength regime. Here, we test expressions (|2|, ^ and (47) against 
the results of our numerical code. 

In Fig. [6] the numerically-determined phase shifts are compared with approximation 
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Figure 2. Absorption cross section of scalar, spinor, EM and gravitational waves up to 
Moj = 0.5. This plot compares the Schwarzschild absorption cross section of massless 
scalar s = 0, spinor s — 1/2, electromagnetic s = 1 and gravitational s = 2 waves. The 
electromagnetic cross section is reproduced from a paper by Crispino, Oliveira, Higuchi and 
Matsas with many thanks to the authors. 



(47) for Muj = 0.05. As expected, there is good agreement. The lower panel shows the 
difference between the numerical and approximate values. It certainly seems plausible that 
the approximation (47) is valid up to 0{e^) ~ (9(0.001) in this case. Note that there is some 



ambiguity in the overall (/-independent) phase factor, which arises from the freedom to choose 



a constant of integration when defining r^, from (27). In these plots, the overall phase term 
has been fixed by matching approximate and numerical phases at / = 2. 

Figure [7] compares the numerically-determined Schwarzschild cross sections at Mu = 0.05 
and Mu = 0.1 with approximation Again, as expected, we find good agreement. Note 
that the helicity-reversing cross section M~'^\g\'^ tends to unity in the backwards direction 

{e = n). 

Figure |8] illustrates the polarizing effect of black hole rotation in the long-wavelength 
regime. The cross sections were calculated by summing the partial wave series up to 



/max = 500 using two iterations of the series reduction technique (5.3) and the asymptotic 



phase shifts (47). It is clear that there is a helicity-rotation coupling which splits the co- 
rotating and counter-rotating helicities. The numerically-determined cross sections are in 
excellent agreement with the asymptotic result ([3]). 

6.2.2. Schwarzschild scattering cross sections. Let us now examine scattering from a non- 
rotating hole in the intermediate regime A ~ r5 {Mu ~ 1). In the following plots, the cross 
section at small scattering angles 9 < 30° is suppressed. At small angles, the series reduction 
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Figure 3. Absorption cross section for a co-rotating wave (to > 0) incident on a fast- 
rotating (a = 0.99M) black hole along the axis of rotation. The top plot shows the 
gravitational absorption cross section CTq as a function of gravitational coupling Muj. The line 
at ctq « 23AttA'P shows the geometric-optics value. The broken curves show the contributions 
CT; from each angular mode, from I = 2 onwards. The bottom plot shows the reflection factor 
for each angular mode. Note the superradiancc (induced emission) exhibited by the I ~ 2, 
m — 2 mode. 



technique (see Sec. 5.3) is sensitive to small numerical errors in the phase shifts at large-Z. 
This sensitivity is due to the divergence on-axis which arises because as / — oo the phase shift 
diverges logarithmically, exp{2iSi) oc Since large-/ partial waves don't 'see' the rotation 

of the hole, rotation does not influence the form of the divergence as 6 ^ 0. At sufficiently 
small angles, approximation ^ is valid. 

Figure |9] shows Schwarzschild scattering cross sections up to Mu = 0.8. As the coupling 
increases, we observe the development of additional structure at large scattering angles. A 
glory halo develops in the backward direction, and at large scattering angles approximation 
([2]) is no longer sufficient. 

Figure 10 shows Schwarzschild scattering cross sections up to Mu = 4. The regular 
oscillations at large angles seen in this figure are predicted by a WKB analysis and semi- 
classical arguments [28l|30]. The angular width of the oscillations is inversely proportional to 
Muj. As expected, the accuracy of the glory scattering approximation ^ improves as Mu 
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Figure 4. Absorption cross section for co-rotating and counter-rotating waves incident on a 
fast-rotating (a* = 0.99M) hole along the axis of rotation. This plot shows that more of the 
counter-rotating wave is absorbed. The middle curve shows the absorption cross section of a 
linearly-polarized incident wave. The line at ctq « 23AttM^ shows the geometric limit [40] . 
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Figure 5. Absorption cross section for a range of black hole rotation rates, a* = 0, 0.5, 0.8, 
0.9, 0.95, 0.99 and 0.999. The left plot shows the cross section for the co-rotating helicity 
(a; > 0) and the right plot shows the counter-rotating helicity (uj < 0). Superradiance 
(ctq < 0) is present only in the co-rotating case. 
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Figure 6. Phase shifts (P = —I) at low coupling Muj — 0.05. The top plots compare the 
numerically-determined phase shifts (defined in Eq. 30 1 with the Iow-Mlu approximation 
(Eq. |47|. The left panel shows the — phase shifts and the right panel shows a* — 0.99. 
It is no surprise to find the phase shifts are only weakly dependent on a* in the low-frequency 
regime, since (47l predicts 5;'^„(a = 0) — 5;„^(a = a*M) ^ 0(a*(Mcj)^). The lower 
panel shows the difference between the numerical and approximate values; the difference 
is consistent with the error 0{M^uj^) in approximation (47 1. 



increases. 

Figure [TT] compares the Schwarzschild scattering cross section of massless waves of spin 
s = (scalar), |s| = 1/2 (spinor/neutrino) [TH] and \s\ = 2 (gravitational). For waves with 
non-zero spin (s > 0), the cross section is zero in the backward direction (note that the 
amplitude g{6) is negligible for large Mu). The scalar wave has a glory maximum in the 
backward direction. The spiral scattering oscillations are of similar width. The first peak of 
the gravitational wave near ^ ~ vr coincides with a peak of the spinor wave. At intermediate 
angles (45° - 120°), the gravitational peaks approximately coincide with the scalar peaks. 

6.2.3. Kerr scattering: long wavelengths Muj < 1. Let us now present results for the 
scattering of waves with a well-defined helicity by a rapidly-rotating (a = 0.99M) black 
hole. The incident wave is travelling along the rotation axis of the hole. Hence, the helicity of 
the initial wave is either aligned {u > 0, 'up', co-rotating, prograde) or anti-aligned {u < 0, 
'down', counter-rotating, retrograde) with the rotation of the hole. We will see that the 
cross section depends on the initial hehcity, and that the scattering process induces a net 
polarization in an unpolarized beam. 
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Figure 7. Schwarzs child scattering cross sections at low couplings. These plots compare 
the numerically-determined scattering cross sections at Mlo = 0.05 and Mlo — 0.1 with the 
approximation Q. Note the logarithmic scale on the y-axis of the lower plot. The thin black 
lines show the helicity-conserving and helicity-reversing contributions. 




Figure 8. Kerr scattering cross sections at low coupling. This plot shows the hclicity- 
preserving |/p and helicity-reversing j^p cross sections at Muj = 0.1 on a logarithmic scale 
for three cases: (i) o = [solid], (ii) a — 0.99M, co-rotating helicity a; > [dashed], and (iii) 
a — 0.99M, counter-rotating helicity u; < [dotted]. The cross sections were numerically 
computed using phase shifts (47 1 and two iterations of the series reduction method (5.3 1. 
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Figure 9. Schwarzschild scattering cross sections (a = 0) at intermediate couplings, 
0.2 < Mu < 0.8. 
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Figure 10. Schwarzschild scattering cross sections (a — 0) at higher couplings 1 < Mto < 4. 
The dotted blue line shows the glory approximation of Eq. [5j 



Scattering and absorption of gravitational plane waves by rotating black holes 



29 



100 



C3 

T3 

"5 
■a 



50 



1 1 


1 1 1 1 

o — U 


1 


cninnr q 


1/2 




gravitational, 


r, O 




















/ ^ 


















1 1 




v\ -'\ » '■- = 1 



45 



60 



75 



90 



105 120 
angle 8 (deg) 



135 150 165 180 



C3 

T3 

O 1 



° 
O 



1 1 



.2 I I I I I I I I I : I 

45 60 75 90 105 120 135 150 165 180 

angle (deg) 



Figure 11. Scattering cross sections for massless waves of spin \s\ — (solid), \s\ ~ 1/2 
(dashed) and \s\ = 2 (dotted) at Mcu — A, a = 0. Note the log scale on the y-axis of the lower 
plot. 



Figure 12 compares the (real part of) the phase shifts of 'up' and 'down' hehcities, for 
a range of angular modes 2 < Z < 60 at Mu = 2. The lower plot shows the proportion of 
each Z-mode that is scattered rather than absorbed. The plot makes it clear that there is a 
critical value l^- below, / < 1^, modes are absorbed; above, / > Ic, modes are reflected. Here, 
Ic depends primarily on Muj, but also on a* and the helicity of the incident wave. Ic is smaller 
for the co-rotating helicity than the counter-rotating helicity. The most significant differences 
in the 'up' and 'down' phase shifts occur close to Ic- At high / ^ Ic, the difference between 
the phase shifts is negligible. In the short-wavelength regime, the critical value Ic may be 
estimated via semi-classical arguments 



Z J UJ z 

where be is the critical impact parameter defined in ([T]). Numerical values for be are listed 
in Table [TJ The value of Ic determines the angular width of the spiral scattering and glory 
oscillations 133). 
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Figure 12. Phase shifts and reflection factors for 2 < I < 60 at Mlo — 2. The upper 
plot compares the phase shifts for co-rotating [Re((5;'^^^g), + symbol] and counter-rotating 
helicities [— R.e((5;'^^^Q), x symbol] for a fast-rotating hole (a = 0.99M) with the Schwarzschild 
phase shifts [* symbol] . The line shows the approximation ( 47 1 (with an arbitrary overall 
phase). The lower plot shows the proportion of each mode which is reflected. 



Figure 13 compares the 'up' (red) and 'down' (blue) scattering cross sections across a 
range of couplings 0.2 < \Muj\ < 0.8. Helicity has a substantial effect on the scattering 
pattern, particularly at large angles > 30°. Superradiance from the Z = 2, m = 2 mode 
enhances the amount of co-rotating flux that is scattered in the backward direction. 

6.2.4- Flux in the Backward Direction 9 = 180°. The backward-scattered flux arises entirely 



from the helicity-reversing amplitude g{6 = vr) (Eq. 16), and hence the backscattered 
component has the opposite helicity to the incident wave. Rapid rotation a ~ 1 enhances 
the effect. A question naturally arises: what is the maximum possible cross section in the 
antipodal direction 9 = 180°? 



Figure 14 shows the cross section in the antipodal direction M~'^\g{TT)\'^ as a function of 
Muj, for a range of rotation rates a* = 0.0 . . . 0.999. We observe an enhancement of ~ 20.3 
times at Moo ^ 0.76 for a = 0.99M, and an enhancement of ~ 34.8 times at Muj ^ 0.945 for 
a = 0.999M. 
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Figure 13. Kerr scattering cross sections (a = 0.99MJ at couplings, 0.2 < Moj < 0.8. The 
dashed red hne shows the cross section for the co-rotating helicity (w > 0) , whereas the dotted 
blue hne shows the cross section for the counter-rotating hehcity (w < 0). 



The transmission factor for the / = 2, m = 2 mode is also shown in Fig. 14 It is clear that 
superradiance plays a role in increasing the back-scattered flux. Nevertheless, the maximum 
of M~^|5f(7r)p does not necessarily coincide with the superradiant peak. Intriguingly, for 
a = 0.999M, the backscattered component has two maxima, at Mu = 0.78 and Mu = 0.945. 
This behaviour can be understood by examining the phase relationship between the / = 2 
mode and the higher modes / > 2. 

Figure 15 illustrates how the helicity-reversed backscattering amplitude g{TT) is 
constructed from a sum over /, as g{9 = n) = 'Ylii=29i- -^o'^ ~ 0.999M, the contribution 
from g2 is largest; nevertheless, the contribution from the higher modes gi^2 is non-negligible, 
and the phase relationship between I = 2 and the higher modes is responsible for the dual- 
peak structure seen in Fig. 14 At Mu = 0.78 and Mu = 0.945 the contributions are 
roughly in-phase (i.e. the arrows in Fig. 15 are roughly aligned), whereas at Mu = 0.88, 



the contributions are roughly out-of-phase (i.e. arrows roughly anti-aligned) . The transition 
from constructive {Mu = 0.78) to destructive {Mu = 0.88) interference and back again 
{Muj = 0.94) occurs because the phase of g2 changes more rapidly with Muj than the phase 



of 5'i>2- This is shown clearly in the right plot of Fig. 15 



6.2.5. Kerr scattering: short wavelengths Muj > 1. Cross sections at higher couplings 
1 < Muj < 4 are shown in Figure 16 Above the critical superradiant frequency, Muc = a/r^, 



the flux scattered in the exact backward direction is negligible. The angular width of the 
spiral scattering oscillations is proportional to 1/lc] hence the 'up' helicity leads to wider- 



Scattering and absorption of gravitational plane waves by rotating black holes 



32 




1 

0.5 


-0.5 
-1 
-1.5 



1 


^^^^^^ — — ''^^^''''^ 


/ / 


1 

/ 






/ / 




1 


1 1 1 




1 



o 
o 
nj 



O 

05 



0.2 



0.4 0.6 
M 0) 



0.8 



Figure 14. Scattering cross section in the backwards direction for prograde wave. The top 
plot shows the flux scatted in the direction opposite to incidence, AI^^^{9 — n) as a, function 
of Mu!, for a co-rotating wave. A range of rotation rates a = 0, . . . 0.999 are plotted. The 
bottom plot shows the transmission factor for the I — 2,m ^ 2 mode. 



angle oscillations than the 'down' helicity. The interaction between BH rotation and wave 
helicity decreases in significance as Muo oo. 

In Figure 17 we examine the glory more closely for the special case Muj = 4.0, a = 0.99M. 
The numerically-determined co-rotating and counter-rotating cross sections are shown as solid 
lines. The broken lines show the semi-classical glory approximations ([5]). The dotted line 
was found using the logarithmic approximation to the scattering angle (13); the dashed line 
using exact solutions to the geodesic equations ([T]). The dot-dash line shows that the glory 
approximation is improved if we include the second-order contribution from rays which orbit 
the hole one-and-a-half times. A significant deficiency of the approximation is that it does 
not take into account the helicity-rotation coupling that splits the co- and counter-rotating 
cases. In future, it would be interesting to model how geodesies are defiected by the effect of 
a spin-rotation coupling. 



6.2.6. The extremal rotation limit, a^, ^ 1. Numerical accuracy becomes harder to obtain 
as a^, ^ 1. The main source of error is in the numerical calculation of phase shifts. A small 
parameter u = (r+ — t_)/2 ^ M 1^2(1 — a^,) appears in the denominator in the near-horizon 
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Figure 15. Back-scattering amplitudes for a ~ 0.999Af . The complex amplitude g(9 = tt) 
arises from a sum over modes, g — X]/=2 5''- T^^*^ left plot shows g (the dotted red arrow), 
and (72 (the long black arrow) and g;>2 (shorter black arrows) in the Argand diagram for 
Muj = 0.72, 0.78, 0.88 and 0.94. The right plot shows g2 (red and blue arrows) and J2i=3 9i 
(green and purple arrows) for a range of couplings Muj — 0.7 . . . 0.98. The phase of 32 
changes more rapidly than the phase of higher modes gi>2, leading to interference effects and 



the double-peak structure in Fig. 14 
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Figure 16. Kerr scattering cross sections (a = 0.99M^ at higher couplings (1 < Muj < 4^. 
The dashed red (dotted blue) line shows the cross section for the co-rotating (counter-rotating) 
helicity with a; > (w < 0) . 
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Figure 17. Glory peaks for AI\u!\ = 4 and a = 0.99M. The solid lines show the partial- 
wave cross sections for the co-rotating {uj > 0) and counter-rotating {u; < 0) helicities. The 



broken lines show the semi-classical approximations derived in Sec. 2.2 The most accurate 
approximation (labelled "2nd order" ) is found by including the additional contribution from 
rays which axe scattered by 6* ~ Stt. 



series expansions. To test the accuracy of the phase shifts, we tried changing the arbitrary 
parameters in the numerical code: the matching radii r^, r^id and Too- We found the change 
in phase to be invariant to at least one part in 10~^ for non-extremal a^,. The error was 
confirmed to be less than the thickness of the lines in all the plots we have presented here. 
However, we found results beyond a^: > 0.9995 to be less reliable. A separate analysis of the 
extremal case a^, = 1 would be of interest. 

7. Discussion and Conclusions 



In this study we have shown that it is possible to compute accurate gravitational-wave cross 



sections using a partial- wave series approach. Although the key formulae ( 15 - 19 ) were derived 



three decades ago by Matzner et al. [TOl [HI [3T] , we believe this is the first time that accurate 
results for the rotating case have been presented. 

Early studies [10], IHl |3T] encountered two obstacles to progress. Firstly, the Teukolsky 
equation has a long-ranged potential for perturbations with spin |s| > 1/2. This leads to 
so-called 'peeling' behaviour in the asymptotic solutions which renders the determination of 



phase shifts numerically difficult. Secondly, the partial wave series ( 15 ) has an infinite number 
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of terms which increase in magnitude with /, and it is divergent at = 0. 

The first obstacle is overcome by transforming the radial equation. Matzner and Ryan [10] 
used Press and Teukolsky's transformation |33]; we used a Sasaki-Nakamura transformation 
[50] . Finding the phase shifts is then relatively straightforward - although of course we benefit 
from the rapid advance in computing power over the last thirty years. 

It is overcoming the second obstacle - approximating an infinite, divergent series - that 
is the key to numerical progress. Regardless of computing power, the series must be truncated 
at some I = /max- Yet, the magnitude of the coefficients in the series (15) increases with /. 
Matzner et al. tried decomposing the phase shifts into a long-range 'Newtonian' contribution 
(that varies logarithmically with /) and a short-range remainder (that decays with /). This 
approach works well for scalar waves [33], but it does not seem to work for higher spins [T9] . 
In truncating the series, Matzner et al. introduce numerical error with an angular frequency 
proportional to l//max- See for example Fig. 14 in [TT]. In this work, we circumvented the 
obstacle by employing a series reduction method to improve the convergence properties of the 
series [BD]- 

The results presented here shed new light on three intriguing scattering phenomena: 
polarization, helicity violation, and glory/spiral diffraction. Below we briefly recap the 
important points. 

The coupling between the rotation of the black hole and the helicity of the incident wave 
means that a rotating black hole has a polarizing effect. In other words, since waves of co- and 
counter-rotating helicities are distinguished by a rotating black hole, a partial polarization is 
induced in an initially unpolarized beam. The polarization is a function of scattering angle. 
In the long-wavelength limit, the partial polarization is given by Q. At shorter wavelengths 
the polarization it is a more complicated function of scattering angle (see e.g. Fig. 16), due 
to diffractive effects. 



The parity-dependence of the phase shifts leads to a scattered component M \g 



(Eq. 16) which has the opposite helicity to the incident flux. The helicity-reversed flux 



is largest in the antipodal direction {9 = 180° = tt). Conversely, the helicity-conserved 
flux |/(^)P is zero in the backwards direction. In the long-wavelength limit we showed that 
M~^|(yf(7r)p = l+4auj+0{uj'^). In the non-rotating case, we flnd |5f(7r)p has a maximum of just 
M~^|(yf(7r)p ~ 1.2 at Mui = 0.3 (Fig. [7|, and it is suppressed at shorter wavelengths (Fig. 10). 
For a co-rotating incident helicity, the back-scattered flux is enhanced by superradiance in the 
/ = 2, m = 2 mode, which occurs for frequencies uj < Uc = a/Mr+. We flnd the backscattered 
flux may be enhanced by an order of magnitude or more as the extremal rotation limit is 
approached (a ^ 1). In particular, for a^, = 0.999 we observe a maximum enhancement of 
M~^|(7(7r)p ~ 35 at Mu ~ 0.945. This contrasts with the scalar-wave (s = 0) case, for which 
superradiance has a negligible effect [33] . 

At short wavelengths, A < r5 {Mu > 1), we observe regular peaks and troughs in 



the scattering cross sections (Figs. 10, 11 and 16). This is a diffraction effect, due to the 



interference of paths scattered through angles 6, 27: — 6,2tt + 6 [3Q\. The angular width of the 
oscillations is (approximately) proportional to the ratio of the inner unstable orbit radius to 
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the incident wavelength. Hence, the angular width is roughly proportional to Mu, but also 
depends on rotation and incident helicity. Again, since co- and counter-rotating helicities 
are distinguished by a rotating black hole, the scattered wave acquires a partial polarization 
which depends on scattering angle. 

Unlike light, gravitational radiation is able to 'see through' the complicated and messy 
astrophysical region surrounding a black hole (e.g. accretion disk, jets, synchrotron radiation, 
etc) to probe the near-horizon geometry. Is there any prospect of detecting the gravitational 
radiation that is scattered by black holes? Since gravitational radiation has yet to be detected 
at all, the short answer is: not yet! The best hope in the short term is that experimentalists 
detect the strongest signals generated by highly dynamic, strong-field interactions, such as the 
final merger of black holes. Dynamic (i.e. time-dependent) perturbations will excite black 
hole 'quasinormal modes' (QNMs) |64|- The frequencies and decay times of QNMs would 
provide a clear signature of the underlying black hole parameters M and a^. 

Diffraction patterns like those shown in Figs. 10, 11_ and 16 may provide an alternative 
signature for future detectors. Let us imagine a supermassive black hole (such as Sag. A*) 
illuminated by a strongly-collimated and long-lasting beam of gravitational radiation. Let 
us assume the incident radiation has a short wavelength A <^ and a wide-band spectrum 
Aujq ^ 1/M. Our angle of observation 6 is effectively fixed, but we are free to vary the 
observation frequency. Changing the observation frequency by 6u! ~ 1/M would bring us 
from an interference maximum to a minimum [21] , and, if the BH is rotating, the polarization 
of the signal would be modified. Hence the characteristic signature of BH scattering would 
be a signal whose intensity and polarization oscillates rapidly over a frequency 6uj. 

Two extensions of this work suggest themselves. Firstly, one could treat waves 
approaching from an arbitary angle of incidence relative to the BH rotation axis. The relevant 
partial wave formulae, given in [31], involve an additional sum over the azimuthal numbers m. 
We would expect to find maximal polarization for waves impinging along the axis of rotation, 
and zero polarization for waves approaching in the equatorial plane. Secondly, it would be 
relatively straightforward to apply the numerical method detailed here to the scattering of 
electromagnetic waves from a rotating hole, since the partial wave formulae are supplied in 
[3T] . and an extension of the Sasaki-Nakamura to |s| = 1 is available in [6^. In addition, it 
should also be straightforward to extend the long-wavelength approximation of section |4] to 
the electromagnetic case. We hope this would shed light on the physical origin of the factor of 
two difference between ^ and the result (38) of Barbieri and Guadagnini [25] for scattering 
from classical rotating matter. We hope to undertake this calculation in the near future. 
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Appendix A. Scattering Angle Approximation 



The left-hand side of ([7| may be expressed in terms of elliptic integrals. For deflection angles 
close to 6 = n, we may expand in powers of a/b, 



bQ = 2K{a/b) - F(cos(^/2), a/6) 
462 + 64 64 ^'^ 



+ 



48 62 



+ . 



(A.l) 
(A.2) 



where K{}i) and /c) are complete and incomplete elliptic integrals of the first kind. 
The right-hand side of ([T]) may be expressed as 

dr 2F{z,k) 



(A3) 



where 



r2 - rp 
rs-ro'' 



and k 



r2 - ri 



rs - ro 



(A.4) 



We are interested in the strong-deflection limit 6 ~ be- In this case, the largest two roots 
r2 and lie close to r^. To derive a logarithmic approximation similar to Darwin's, let 



r2 
k 



r,{l-5), r3 = re(l + 5). 



-5, z 



where 



6' 



Tc-ro (rc-ri)(r, 
(2r2 - 4Mrc + 2a2) 62 /6 - 6, 



(6r2 + 2a2-62) r2 
Let us now make use of the asymptotic expansion 

F{z,k) In 



+ Oi{b-b,y). 



(A5) 
(A6) 

(A.7) 



1 + (1 + ^'2^2/^/2)1/2 

where z' = \J\ — z'^ and k' = y/1 — k"^. Then the right-hand side of ([7| is approximately 

In 



dr 



(l+^/l4^)4^2g2 
64(rc— ro)2 



R 2^J{rc - ro)(rc -ri) 



(A9) 



where a = {k'/z')^ = (ri — ro)/(rc — ri). Combining (A.2) and (A. 9) we reach the conclusion 
that the scattering angle can be expressed in the logarithmic form (13), where 

6. 



C{a) 



and 
1 



(r,-i?o)(rc-i?i) 
1 + v^TTa) (2r2 - 4Mre + 2a2) Mbc 



D{a) QA{r,-Roy (6r2 + 2a2 - 62) 



exp 



TT 



C{a) 



(A.IO) 



[1 + aV262 + 9a7326^] (A. 11) 
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This section describes a stable method for accurate computation of the spherical harmonics of 
spin weight s = —2 using recurrence relations [51]. The zero-spin- weight spherical harmonics 
are written 



^Yi^{9) = A{l,m)PiUcos9), 



where 



A{l,m) 



21 + 1 {l-m)\ 



An 



(B.l) 



(B.2) 



(/ + m)\ ' 

The associated Legendre polynomials Pim{x) can be computed via the recurrence relation 

(/ - m) = x{2l- l)P,_i,„ - (/ + m - l)Pi-2,m (B.3) 

with initial values 

Pmrnix) = (-l)"(2m - 1)!! (1 - X^)™/^ 

(x) = x{2m + l)Pmm. (B.4) 
The operator 6 lowers the spin weight of a function [66], as follows 



SsYimiO) = -(sin( 



d m 
+ 



36 sin 6 

= -[{l + s){l-s + (s-i)YUe). (B.5) 

Hence the s = —2 functions are found by acting twice with a spin-lowering operator 6 on the 
zero-spin- weight spherical harmonics o^/m- This yields 



A{l,m) 



v/(/ -!)/(/ + !)(/ + 2) 



:i-x' 



I — 
fix" 



- 2m— + 
ax 



d rn^ — 2mx 



x^ 



Pmix) 



(B.6) 



To compute these harmonics, we need recurrence relations for dPim/dx and (1 — x'^)d'^Pim/dx. 

dPl-2,r 



These are straightforward to derive from (B.3) and (B.4): 
dPim 



dx 
dP 



1 

/ — m 



{21 - 1) P,_i,^ + X 



dPl-l^r. 

dx 



(/ + m- 1)- 



dx 



dx 
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mx(-l)'"(2m - 1)!!(1 - x^Y"'-^^/^ 
dP 



dx 



(2m 
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d-'P 



Im 



dx^ 



'1 — X 



I — m 
2 



(2/- 1 



dx 
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2 ^ \-x- 



(l + m-l) 



d'Pi- 



2,m 
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m 
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The coefficients a'^ appearing in the MST formalism (46-45) are 



V 

a„2 
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With the Cordon- Short ley phase convention, the Clebsch-Gordan coefficients appearing 



in (54-59) are 
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The expansion coefficients {d^25 '^-15 '^-i? '^o^^ '^+15 '^+15 '^+2} appearing in (53) are 
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The function F(r) appearing in (102) is given by 



where 



and 



li^) = Co + Ci/r + C2/r^ + Cs/r^ + C4/r^ (D.2) 



Co = -l2iuM + A(A + 2) - I2auj{aw - m), 
ci = 8m [Sau — \{auj — m)] , 

C2 = -2AiaM{auj - m) + 12a^ [l - 2{auj - mf] , 
C3 = 2Aia^{auj - m) - 24Ma^ 

C4 = 12a^ (D.3) 
The function U{r) appearing in ( |102 ) is 



(r2 + a2)2 ^ ^ ' ^ r^ + a'' 



Uir) = + GirY + - F(r)G(r) (D.4) 



where 
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;7i(r) = nr) + ^ 



2 
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The function K(r) was defined in Eq. (24), and V{r) is the Teukolsky potential 
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